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摘 要 : 浑 善 达 克 沙 地 作为 中 国 北 方 沙漠 化 最 严重 的 地 区 之 一 ,已 被 列 入 全 国 重点 生态 环境 建设 
地 区 ,其 沙漠 化 发 展 过程 及 其 影响 因素 研究 具有 重要 意义 。 利 用 1985 一 2017 年 浑 善 达 克 沙 地 5 期 


遥感 影像 数据 ,通过 eCognition 9.0 软件 ,采用 本 


向 对 象 计算 机 自动 分 类 法 ,获取 沙漠 化 指数 ,进行 


浑 善 达 克 沙 地 沙 汉化 时 空 演变 特征 研究 ,并 利用 主 成 分 分 析 法 ,进行 了 沙 汉化 驱动 机 制 研 究 。 结 


果 表 明 : 近 30a, 浑 善 达 克 沙 地 沙漠 化 土地 
轻 度 沙漠 化 土地 


积 以 2002 年 为 转折 点 呈现 先 增加 后 减少 的 趋势 ,其 中 
可 先 减 少 后 增加 ,由 42.19% 减 少 到 20.9% 后 又 增加 到 32.1% , 主要 分 布 于 东部 地 


区 ;中 度 沙 漠 化 有 逐渐 增加 的 趋势 ,由 20.6% 增 加 到 39.0% ,分 布 于 中 部 地 区 ;重度 沙漠 化 土地 由 
30.3% J > Bl] 22.6% , 而 极度 沙漠 化 先 增 加 后 减少 ,由 7.0% 增 加 到 16.0% 后 又 减少 到 6.3% ,分 布 于 


西部 及 以 南 逐 渐 往 西北 方向 移动 。 沙 汉化 程度 空间 上 形成 由 西向 东 逐 渐 减 轻 的 趋势 。 浑 善 达 克 
沙 地 沙漠 化 发 展 过 程 受 年 末 总 人 口 、 国 内 生产 总 值 . 三 次 产业 、 载 畜 量 、 耕 地 


只 等 人 为 因素 影响 


较 大 ,自然 因素 起 到 辅助 作用 。 因 此 , 沙 汉化 治 到 


过程 中 ,有 效 实施 退耕 还 林 还 草 和 禁 牧 \、 休 牧 等 


政策 ,合理 开展 治 沙 防 沙 工 程 ,进一步 改善 浑 善 达 克 沙 地 沙漠 化 ,推动 区 域 生态 环境 的 健康 发 展 。 
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文章 编号 : 


1977 年 联合 国 召 开 沙 漠 化 会 议 以 来 ,沙漠 化 作 
为 一 个 生态 环境 问题 被 世界 所 公认 "。 长 期 以 来 ， 
全 世界 110 多 个 国家 受 沙漠 化 危害 ,而 中 国 是 世界 
上 上 受 沙漠 化 危害 最 严重 的 国家 之 一 ,分 布 广 .面积 
大 是 其 最 主要 的 特点 ”。 沙 漠 化 对 中 国 北方 干旱 E 
干旱 地 区 生态 环境 和 社会 经 济 的 影响 尤为 重要 **。 
沙漠 化 ( 沙 质 荒漠 化 ) 作 为 极端 王 旱 .干旱 TEA 
部 分 羊 湿 润 地 区 土地 退化 过 程 的 典型 特征 之 一 ， 
对 所 在 区 域 生态 .水文 和 农业 等 方面 均 造 成 严重 危 
害 , 如 沙尘暴 发 生 .土壤 风蚀 草场 退化 .耕地 生产 
力 下 降 等 ,从 而 严重 影响 区 域 社会 经 济 的 发 展 ”。 
因此 ,一 直 以 来 沙漠 化 问题 不 仅 是 社会 各 界 广 泛 关 
注 的 焦点 ,其 变化 .成因 防治 等 方面 ,也 是 地 理学 、 
生态 学 等 学 科研 究 的 热点 之 一 。 


收 稿 日 期 : 2020-05-19; 修订 日 期 : 2020-07-24 


20 世 纪 50 年 代 , 中 国 科 学 院 初 步 开 展 了 有 关 沙 
漠 化 的 研究 工作 ”。1978 年 , 朱 震 达 等 ”首次 提出 
“沙漠 化 ”的 概念 ,并 对 其 定义 和 形成 机 理 做 了 系统 
的 阐释 。 到 了 20 世纪 80 年代 ,中 国学 者 们 逐步 开展 
了 沙 演化 相关 研究 工作 ,主要 围绕 沙漠 化 的 成 因 ”、 
空间 格局 、 演 变 过 程 ” 防治 措施 “和 脆弱 性 
评估 “等 方面 ,其 中 沙漠 化 动态 监测 是 其 他 研究 的 
基础 。 随 着 遥感 技术 的 发 展 ,其 对 地 表 过 程 的 动态 
监测 方面 起 到 了 很 重要 的 作用 ,从 而 在 沙漠 化 研究 
方面 也 成 为 主要 的 监测 手段 ”"。 通 过 遥感 手段 ,对 
中 国 典型 的 沙漠 化 地 区 ,如 内 蒙古 阿拉 善 高 原 时 
尔 多 斯 高 原 “、 浑 善 达 克 沙 地 科尔沁 沙 地 “以 
及 青藏 高 原 ”和 青海 共和 盆地 ”等 沙漠 化 地 区 , 进 
行动 态 监测 的 基础 上 研究 其 时 空格 局 ,从 而 进一步 
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研究 其 形成 机 制 。 而 沙漠 化 驱动 因子 是 分 析 其 演 
变 和 防治 的 前 提 , 通 和 常 沙漠 化 的 影响 因素 主要 有 自 
然 因 素 和 人 为 因素 两 大 类 ,采用 数理 统计 方法 研 
究 驱 动因 子 与 沙漠 化 之 间 的 关系 。 常 学 礼 等 ” 采 
用 累加 效应 分 析 法 研究 自然 因素 和 人 为 因素 对 科 
尔 沁 沙 地 沙漠 化 的 影响 ,进一步 利用 3S 技 术 分 析 了 
典型 地 区 人 工 造 林 对 沙漠 化 过 程 的 影响 当 ; 赵 哈 林 
等 中 王 静 妆 等 3 利用 数理 统计 等 方法 分 别 研究 了 
气候 变化 对 北方 农 牧 交错 区 和 科尔沁 沙 地 土地 沙 
漠 化 的 影响 ; 阿 如 旱 等 利用 Logistic 回归 模型 分 析 
了 内 蒙古 多 伦 县 土地 沙漠 化 驱动 因子 。 而 已 有 的 研 
究 , 主 要 进行 某 一 时 间 点 的 驱动 因子 和 沙漠 化 两 者 
之 间 的 关系 研究 ,并 采用 沙漠 化 的 面积 ”指数 字 3 
作为 沙漠 化 程度 的 评估 指标 ,与 已 有 研究 相 比 ,本 
文 数据 有 所 更 新 .研究 跨度 有 所 增加 (1985 一 2017 
年 ), 且 进行 多 期 驱动 因子 与 沙漠 化 关系 研究 。 

浑 善 达 克 沙 地 处 于 中 国 北方 典型 草原 区 向 旱 
作 农 业 区 过 渡 地 带 , 是 中 国 北方 荒漠 化 最 严重 的 地 
区 之 一 ””。 作 为 中 国 沙尘暴 主要 沙 尘 源 地 之 一 ， 
家 对 浑 善 达 克 沙 地 实施 了 退耕 还 林 还 草 . 禁 牧 休 
牧 等 治理 沙 尘 源 和 生态 环境 恢复 的 一 系列 政策 ,使 
该 地 区 沙漠 化 发 展 势头 得 到 过 制 ,区域 生态 环境 有 
所 改善 。 但 是 , 浑 善 达 克 沙 地 生态 环境 脆弱 ,加 
之 气候 变化 较 大 ,不 合理 的 土地 利用 方式 和 过 度 开 
暴 、 放 牧 等 现象 并 没有 得 到 彻底 改善 ,导致 该 地 生 
态 环境 仍然 面临 问题 ,从 而 直接 影响 了 社会 经 济 的 
可 持续 发 展 ? 。 因 此 ,本文 以 浑 善 达 克 沙 地 为 研究 
区 ,利用 1985、1992、2002、2011 年 和 2017 年 等 5 期 
的 Landsat TM/OLI 遥感 影像 数据 为 数据 源 ,通过 
eCognition 9.0 软件, 采用 面向 对 象 计算 机 自动 分 类 
法 ,以 沙漠 化 差 值 指数 (DDI) 作 为 衡量 沙漠 化 程度 
的 指标 ,获取 研究 区 沙漠 化 遥感 监测 数据 ,并 利用 
主 成 分 分 析 方 法 ,进行 沙漠 化 变化 过 程 与 自然 .人 
为 因素 之 间 的 关系 研究 ,进而 确定 了 近 30 a 来 浑 善 达 
克 沙 地 沙漠 化 的 时 空 演 化 过 程 和 驱动 机 制 ,为 相关 部 
门 进行 沙漠 化 防治 和 管理 等 工作 提供 科学 借鉴 。 


1 研究 区 概况 
浑 善 达 克 沙 地 地 处 内 蒙古 自治 区 锡林郭勒 大 


草原 南端 ,东部 从 大 兴安 岭 的 南 段 西 坎 ,西部 延伸 
至 集 二 铁路 沿线 以 西 ,东西 宽 约 473 km ,南北 宽 50~ 


100 km, ,地 理 位 置 为 111?27'34.2"~117?10'46.9" E, 
41°10'10.5"~42°58'30.7" N ,总 面积 为 41639.7 km’; 
包括 隶属 与 锡林郭勒 盟 的 苏 尼 特 右 旗 、 共 尼 特 左 
旗 、 阿 巴 嘎 旗 、 镶 黄旗 TEE ATE 锡林浩特 市 .多 伦 
县 和 正 蓝 旗 等 8 个 旗 县 和 赤峰 市 内 的 克 什 克 腾 旗 1 
个 旗 县 (市 ) (图 1); 属 于 中 温带 半 干 旱 .干旱 大 陆 
性 季风 气候 区 ” , 据 1985—2017 年 12 个 气象 站 观 
测 数据 ,年 平均 气温 1.5~5.0 %C ,年 降水 量 200~400 
mm, 自 西 北向 东南 逐渐 增多 ,主要 集中 在 7 一 9 月 ， 
年 均 风 速 2.7~3.5 ms ,年 最 大 风速 6.5~8.0 ms ，， 
夏季 炎热 少雨 ,冬季 寒冷 多 风 ; 地 带 性 土壤 为 栗 钙 
土 , 非 地 带 性 土壤 主要 有 风沙 土 . 盐 碱土 和 草 馈 土 ; 
植被 以 大 针 茅 (Stipa grandis) . 56 ETH (Stipa krylo- 
vii) $Æ (Leymus chinensis ) YKFE (Agropyron crista- 
tum) 地 带 性 植被 和 榆树 (Ulmus pumila) 、 锦 鸡 儿 
(Caragana sinica) 、 树 柳 (Tamarix chinensis) 、 黄 柳 
(Salix gordejevii) , 2 E WE jes (Artemisia halodendron ) 
等 沙 地 植被 为 主 。 


2 数据 来 源 和 研究 方法 


2.1 数据 来 源 和 预 处 理 

2.1.1 RAB RA ”沙漠 化 遥感 监测 数据 来 源 于 
Landsat5 TM/Landsat8 OILT 传 感 咒 数据 ,由 美国 地 质 
调查 局 (United States Geological Survey , https://glovis. 
usgs.gov/) 下载 获取 ,根据 最 早 数据 的 可 获得 性 可 把 
溯 到 1982 年 ,同时 考虑 约 10 a 的 时 间 间 隔 ,分 别 选 
取 1985、1992、2002、2011 年 4 期 的 Landsat 5 TM 数 
据 ( 多 光谱 波段 分 辨 率 30 m) 和 2017 年 的 Landsat 8 
0OLI 数 据 ( 多 光谱 波段 分 辨 率 30m ,全 色 波 段 分 辩 率 
15 m)。 遥 感 影像 主要 选取 植被 生长 季 (7 一 9 H), 
且 云 量 覆盖 小 于 10% 的 影像 ,并 补充 了 少数 非 生长 
季 和 其 他 可 以 反映 不 同类 型 地 物 差异 的 影像 数 
据 。 和 覆盖 整个 浑 善 达 克 沙 地 需要 6 景 数据 ,获取 
1985—2017 年 5 个 时 期 Landsat TM/OLI 数据 共 
30 景 。 

使 用 ArcGIS 10.3 和 ENVI 5.3 软件 对 数据 进行 
几何 校正 、 辐 射 校正 .大气 校正 .投影 转换 和 馈 般 等 
预 处 理 , 并 利用 浑 善 达 克 沙 地 矢量 边界 进行 数据 

沙漠 化 遥感 监测 分 类 方法 采用 面向 对 象 计算 
机 自动 分 类 方法 ” ,并 以 沙漠 化 差 值 指数 (DDU) 作 
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苏 尼 特 右 旗 


图 例 
e 气象 站 点 

O 浑 善 达 克 沙 地 界线 
O 浑 善 达 克 沙 地 相关 旗 县 界 


0 25km 
本 = 


All 浑 善 达 克 沙 地 示意 图 
Fig.1 Hunshandake Sandy Land range map 


为 衡量 沙漠 化 程度 的 指标 ,获取 研究 区 沙漠 化 遥感 
监测 数据 ,在 此 基础 上 结合 野外 调查 和 高 分 一 号 
(GF1) 数 据 进行 数据 结果 验证 ,精度 达到 90% 以 上 。 
2.1.2 驱动 因子 数据 1985 一 2017 年 浑 善 达 克 沙 地 
气象 数据 从 国家 气象 信息 中 心 获取 ,包括 年 平均 气 
温 、 年 降水 量 、 年 平均 风速 和 年 最 大 风速 等 指标 ; 
1986 一 2017 年 的 人 口 数 据 eh at RT A IN 
生产 总 值 和 三 次 产业 等 数据 主要 由 内 蒙古 历年 统 
计 年 鉴 获取 。 

2.1.3 沙 江 化 等 级 划分 根据 Liu 等 “对 沙漠 化 指 
数 的 判别 标准 和 分 类 方法 ,将 1985 一 2017 年 不 同时 
期 浑 善 达 克 沙 地 的 沙漠 化 指数 划分 成 非 沙漠 化 , 轻 
度 .中 度 .重度 .极度 等 $ 个 等 级 ,如 表 1 所 示 。 


表 1 沙漠 化 等 级 划分 表 


Tab. 1 Desertification classification table 


类 别 SDI 值 
极度 沙漠 化 土地 SDI<0.3730 
重度 沙漠 化 土地 0.3730<SDI< 0.4965 
中 度 沙 漠 化 土地 0.4965<SDI< 0.6414 
轻 度 沙漠 化 土地 0.6414<SDI< 0.8024 


非 沙漠 化 土地 
注 :SDI 为 沙漠 化 指数 。 


0.8024<SDI< 1.0977 


2.2 研究 方法 
通过 面向 对 象 分 类 方法 ,提取 沙漠 化 相关 数 
据 ,该 方法 除了 考虑 影像 的 波段 纹理 等 信息 以 外 ， 
还 有 陆地 和 水 体 指数 . 归 一 化 植被 指数 、Albedo 等 
表征 特定 地 物 的 指数 。 
2.2.1 陆地 和 水 体 指 数 (Land and water index, 
LWI) 
LWI= Infra - RED/(GREEN +0.0001)x 100 (1) 
式 中 :Infra-RED GREEN 分别 为 红外 、 绿 光波 段 ,分 
别 对 应 Landsat TM 图 像 的 第 5、2 波 段 和 Landsat OLI 
图 像 的 第 6.3 波 段 。 
2.2.2 沙 江 化 差 值 指数 (Desertification difference 
index, DDI) NDVI 与 植被 覆盖 度 成 正比 ,而 Albe- 
do 与 植被 覆盖 度 成 反比 , 即 Albhedo 和 NDVI 两 者 间 
呈 负 相关 关系 ,说 明 NDVI 值 的 减少 使 沙漠 化 程度 
增加 ,而 Albedo 值 的 增加 使 沙漠 化 程度 也 随 着 增 
加 , 即 沙漠 化 程度 可 以 由 Albedo NDVI 的 一 元 线性 
回归 表达 式 反映 。DDI 可 以 构建 为 
DDI = 2.234 x NDVI-Albedo (2) 
式 中 :NDVI 由 公式 (3) 计 算 ,Albedo 由 公式 (4)、(5) 
计算 。 


202108.00030v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


同 丽 嘎 等 : 近 30 a 浑 善 达 克 沙 地 沙漠 化 时 空 演变 特征 及 驱动 机 制 研究 


2.2.3 归 一 化 植被 指数 (Normalized difference veg- 
etation index, NDVI) NDVI 主 要 应 用 于 监测 植被 
生长 状态 .植被 覆盖 度 和 消除 部 分 辐射 误差 等 ,ND- 
VI 计算 公式 为 : 
NDVI =(NIR-RED)/(NIR + RED) (3) 
式 中 :NIR 和 RED 分 别 为 近 红 外 、 红 光波 段 , 分 别 对 
Di Landsat TM 图 像 的 第 4. 3 波段 和 Landsat OLI 图 像 
的 第 5.4 波 段 。 
2.2.4 Albedo Landsat TMOLI 反 照 率 的 转化 使 用 
的 是 Liang °°" RSE ETT ,公式 如 下 : 
Albedo(TM) = 0.356 x b1 + 0.130 x b3 +0.373 x 
b4 + 0.085 x b5 + 0.072 x b7 — 0.0018 
Albedo(OLI) = 0.356 x b2 + 0.130 x b4 + 0.373 x 
b5 + 0.085 x b6 + 0.072 x b7 — 0.0018 
式 中 :bl、.b2、b3、b4、b5 、.b6、b7 分 别 对 应 Landsat TM/ 
OLI 的 1-7 波 段 。 
2.2.5 沙漠 化 驱动 因素 分 析 沙漠 化 受 自 然 和 人 为 
双重 因素 的 影响 ,是 多 个 因素 综合 作用 的 过 程 。 自 
然 因素 选取 浑 善 达 克 沙 地 范围 内 旗 县 和 周边 共 12 
个 气象 站 点 1985—2017 年 的 年 降水 量 、 年 平均 气 
温 、 年 平均 风速 .年 最 大 风速 等 气象 指标 ,采用 Arc- 
GIS 10.3 克 里 金 插值 法 对 12 个 站 点 的 数据 进行 差 
值 ,并 通过 统计 计算 ,获取 浑 善 达 克 沙 地 其 他 旗 县 
的 气象 数据 。 人 为 因素 选取 1986—2017 年 的 耕地 
面积 . 载 畜 量 及 人 口 数 量 国内 生产 总 值 和 三 次 产 
业 数 据 等 。 由 于 自然 因素 对 沙漠 化 的 影响 是 个 漫 
长 的 过 程 汪 ,因此 ,针对 自然 因素 进行 了 滑动 平均 
处 理 (1~10 a) ,获取 与 DDI 关 联 度 最 高 的 平均 值 结 
果 。 在 此 基础 上 ,文章 采用 主 成 分 分 析 法 进行 浑 善 
达 克 沙 地 沙漠 化 驱动 机 制 分 析 。 该 方法 不 仅 能 简 
化 多 个 变量 ,还 能 量化 沙漠 化 驱动 因素 中 多 个 因素 
的 贡献 率 ,并 揭示 沙漠 化 过 程 的 主要 驱动 力 “。 文 
章 中 选取 的 主要 自然 和 人 为 影响 因素 为 年 降水 量 
(X1) 年 平均 气温 (X2)、 年 平均 风速 (X3)\ 年 最 大 风 
速 (X4) AEA AC (XS) 、 耕 地 面积 (X6) RAE 
(X7) 国内 生产 总 值 (X8)、 第 一 产业 (X9)、 第 二 产业 
(X10) 和 第 三 产业 (X11)。 
2.3 数据 精度 验证 
文章 采用 高 分 辨 率 数 据 (CF1) 以 及 野外 采样 点 
(30 个 ) 来 评价 浑 善 达 克 沙 地 沙漠 化 遥感 监测 数据 
的 质量 。 以 2017 年 的 数据 为 例 ,根据 野外 采样 点 的 
分 布 情况 ,计算 高 分 一 号 (CCGF1) 影 像 的 NDVI 值 , 然 


(4) 


(5) 


后 转换 成 植被 盖 度 ,根据 沙漠 化 指标 分 级 标准 评价 
各 个 采样 点 土地 的 沙漠 化 程度 。 在 Landsat 影 像 范 
围 选取 各 级 别 沙漠 化 验证 点 ,与 本 数据 集 所 得 到 的 
结果 进行 比较 ,并 计算 Kappa 系 数 。 经 验证 ,2 种 分 
辩 率 的 数据 评价 结果 基本 一 致 ,总 体 精度 达到 90% 
以 上 ,反映 出 本 文 的 数据 在 区 域 尺 度 上 评价 浑 善 达 
克 沙 地 沙漠 化 具有 可 行 性 。 


3 结果 与 分 析 


3.1 浑 善 达 克 沙 地 沙漠 化 时 空 演化 特征 
3.1.1 沙 湛 化 土地 面积 时 室 变 化 分 析 近 30a 浑 善 
达 克 沙 地 沙漠 化 土地 以 2002 年 为 转折 点 总 体 上 呈 
现 出 先 增加 后 减少 的 变化 趋势 ,其 中 轻 度 沙漠 化 土 
地 面积 先 减少 后 增加 ,中 度 沙 漠 化 土地 面积 呈现 逐 
渐 增 加 的 趋势 ,重度 和 极度 沙漠 化 土地 面积 呈现 先 
增加 后 减少 ,与 沙漠 化 总 体 变化 趋势 一 致 (图 2)。 
1985—2002 年间 , 浑 善 达 克 沙 地 沙漠 化 土地 面 
积 从 36940.0 km? 增 长 到 38259.7 km ,其 中 轻 度 沙漠 
化 面积 有 明显 减少 趋势 ,从 占 总 沙漠 化 土地 面积 的 
42.1% Ik FI] 20.9% ,而 中 度 和 极度 沙漠 化 土地 面积 
增加 显著 ,分 别 从 20.6% 、7.0% 增 加 到 33.5% 、 
16.0% , 重度 沙漠 化 土地 面积 有 所 增加 ,从 11195.3 
km 增加 到 11350.9 km ,但 所 占 总 沙漠 化 土地 面积 
比例 反而 有 所 减少 ,从 30.3% 减 少 到 29.7% ;2002 一 
2017 年 间 , 沙 漠 化 土地 面积 从 38259.7 km 减少 到 
35968.2 km ,其 中 重度 .极度 沙漠 化 土地 面积 明显 
减少 ,分 别 从 占 29.7% 、16.0% 减 少 到 22.6% .6.3% , 


ce lM — ’‘’ 
1.4} 
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图 2 1985 一 2017 年 浑 善 达 克 沙 地 沙漠 化 土地 面积 变化 
Fig.2 Changes of desertification land areas of Hunshandake 
Sand Land from1985 to 2017 
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而 重度 .极度 沙漠 化 转 为 中 度 和 轻 度 沙漠 化 ,因此 
轻 度 .中 度 沙 漠 化 面积 增长 显著 ,分 别 增长 3542.7 
km’, 1224.6 km ,所 占 比例 分 别 增长 到 32.19% 、 
39.0%. 

浑 善 达 克 沙 地 沙漠 化 程度 空间 分 布 呈 现 由 西 
向 东 逐 渐 减 轻 的 趋势 ,重度 沙漠 化 土地 主要 分 布 于 
西部 ( 苏 尼 特 右 旗 中 东部 . 苏 尼 特 左 旗 西 南部 TE 
白旗 西部 地 区 和 灸 黄旗 北部 小 部 分 地 区 ) ,极度 沙 
漠 化 土地 从 西部 以 南 逐 渐 往 西北 方向 移动 ( 苏 尼 特 
右 旗 中 至 东部 ) ,中 度 沙 漠 化 土地 主要 分 布 于 中 部 
地 区 ( 正 久 白旗 中 北部 ` 苏 尼 特 左 旗 东 南部 、 阿 巴 嘎 
旗 南 部 和 正 蓝 族 西 部 地 区 ) , 轻 度 沙漠 化 土地 主要 
分 布 于 东部 地 区 ( 正 蓝 旗 中 至 东部 大 部 分 地 区 、 锡 
林 浩 特 市 南部 小 部 分 地 区 、 克 什 克 腾 旗 西部 和 多 伦 
县 西北 部 等 地 区 ) ,而 非 沙 漠 化 土地 分 布 于 东部 和 东 
北 地 区 ( 科 什 克 腾 旗 中 部 和 多 伦 县 东部 地 区 )( 图 3)。 


(a) 1985 年 


图 例 


] 轻 度 沙漠 化 土地 m 非 沙 漠 化 土地 
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(c) 2002 年 
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图 例 
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(e) 2017 年 


图 例 


瑟 轻 度 沙漠 化 土地 国 非 沙漠 化 土地 


3.1.2 沙 党 化 土地 转移 变化 分 析 ”根据 浑 善 达 克 沙 
地 沙漠 化 转移 矩阵 ( 表 2) ,1985 一 2017 年 间 , 以 2002 
年 为 转折 点 浑 善 达 克 沙 地 沙漠 化 面积 呈现 先 增加 
后 减少 的 趋势 ,其 中 极度 .重度 沙漠 化 分 别 逆向 转 
化 为 重度 .中 度 沙漠 化 面积 明显 要 大 于 正 向 的 转 
化 。 其 中 1985 一 2002 年 间 , 轻 .中 .重度 沙漠 化 分 别 
转化 为 中 、 重 和 极度 沙漠 化 面积 达到 9031.1 km’, 
5117.1 km’ 、5065.6 km’ , 非 沙 漠 化 土地 面积 减少 达到 
3469.4 km。 显 然 , 浑 善 达 克 沙 地 低级 别 沙漠 化 往 
高 级 别 逐 级 转化 现象 非常 显著 ,沙漠 化 面积 扩大 ， 
总 体 有 明显 加 重 趋 热 ;2002 一 2017 年 , 重 、 极 度 沙 漠 
化 都 有 所 减轻 ,重度 变 为 中 度 沙 漠 化 面积 达到 
8396.3 km ,极度 转化 为 重度 沙漠 化 土地 面积 为 
4514 km , 轻 度 沙漠 化 土地 面积 增加 2124.1 km’, JE 
沙漠 化 土地 面积 增加 1002.1 km ,该 时 期 浑 善 达 克 
沙 地 沙漠 化 程度 有 明显 好 转 趋势 。 


(b) 1992 年 
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(d) 2011 年 
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T 轻 度 沙漠 化 土地 
Bw 中 度 沙漠 化 土地 国 居民 
国 重度 沙漠 化 土地 æ 水 体 0 50km 


ul 极度 沙漠 化 土地 


国 | 中 度 沙漠 化 土地 国 居民 地 


m 重度 沙漠 化 土地 me 水 体 
极度 沙漠 化 土地 


0 50 km 
Ly 


图 3 1985 一 2017 年 浑 善 达 克 沙 地 沙漠 化 空间 格局 图 
Fig. 3 Spatial pattern of desertification in Hunshandake Sandy Land from 1985 to 2017 
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同 丽 嘎 等 : 近 30 a 浑 善 达 克 沙 地 沙漠 化 时 空 演变 特征 


表 2 1985—2017 年 浑 善 达 克 沙 地 沙漠 化 转移 矩阵 


Tab. 2 Transition matrix of desertification in Hunshandake Sandy Land from 1985 to 2017 /km 
类 型 轻 度 沙漠 化 土地 中 度 沙 汉化 土地 ”重度 沙漠 化 土地 ”极度 沙漠 化 土地 非 沙漠 化 土地 总 计 
轻 度 沙漠 化 土地 9282.6 3802.8 233.0 2.2 2239.1 15559.7 
中 度 沙 漠 化 土地 1056.9 5439.5 1002.3 40.1 55.4 7594.2 
重度 沙漠 化 土地 273.4 4341.8 5342.4 1201.7 35.9 11195.2 
极度 沙漠 化 土地 8.3 305.5 1397.5 849.5 30.1 2590.9 
非 沙漠 化 土地 919.4 142.8 136.1 190.4 3311.4 4700.1 
总 计 11540.6 14032.4 8111.3 2283.9 5671.9 41640.1 


近 30 a Æ , 浑 善 达 克 沙 地 沙漠 化 面积 的 变化 主 
要 与 内 蒙古 总 体 荒漠 化 治理 政策 相关 。 从 20 世纪 
50 年 代 开 始 ,内 蒙古 党 委 ,政府 发 布 相关 荒漠 化 防 
治政 策 文件 ,主要 内 容 为 保护 森林 、 防 止 滥 砍 盗 伐 、 
禁止 盲目 开荒 等 。20 世 纪 90 年 代 开 始 ,荒漠 化 防治 
工作 真正 进入 全 面 实施 阶段 ,结合 国家 生态 环境 建 
设 工程 ,从 农业 、 林 业 .畜牧 业 以 及 土地 资源 等 多 方 
面 进 行 荒漠 化 综合 防治 工作 。 因 此 ,从 2002 年 开 

台 , 浑 善 达 克 沙 地 沙漠 化 程度 开始 有 所 好 转 , 且 沙 
漠 化 面积 逐渐 有 减少 的 趋势 。 
3.2 浑 善 达 克 沙 地 沙漠 化 成 因 
3.2.1 沙 党 化 影响 因子 分 析 在 影响 沙漠 化 的 自然 
因子 中 ,气温 、 降 水 量 、 风速 等 气候 因子 与 沙漠 化 之 
间 的 关系 较 密 切口。 近 30 a, 浑 善 达 克 沙 地 年 平均 
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C 年 平均 气温 一 年 平均 风速 
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耕地 面积 /108 hm, 载 畜 量 /108 只 ， 


气温 和 降水 量变 化 波动 较 大 ,总 体 上 呈现 增长 趋 
势 , 其 中 年 平均 气温 为 (3.5+0.8) % ,从 1997 年 开始 
有 所 增长 ,年 降水 量 为 (291.7+51.5) mm, 在 1992 一 
1998 年 和 2012 一 2017 年 2 个 时 期 增长 趋势 明显 。 
而 浑 善 达 克 沙 地 年 平均 风速 为 (3.2+0.2) ms ,年 最 
大 风速 为 (7.3+0.5) m.s"… ,基本 处 于 平稳 状态 (图 4)。 

在 影响 沙漠 化 的 人 为 因素 中 ,主要 选取 了 年 未 
总 人 口 耕地 面积 、 载 畜 量 及 国民 经 济 相关 指标 ,这 
些 指标 作为 沙漠 化 发 生 和 发 展 的 主要 内 容 。 近 30a 
来 , 浑 善 达 克 沙 地 耕地 面积 呈现 先 增 长 (1988 一 
2000 年 ) 后 下 降 (2000 一 2005 年 ) 最 后 趋 于 平稳 变化 
的 趋势 。 而 载 冀 量变 化 波动 较 大 ,总 体 呈 现 先 增长 
后 急剧 下 降 (1999 一 2001 年 ,由 135.5x10!: 只 下 降 到 
87.8x 10° R ) 并 逐渐 趋 于 稳定 状态 。 该 变化 与 1980 
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图 4 1985 一 2017 年 浑 善 达 克 沙 地 沙漠 化 驱动 因子 
Fig. 4 Driving factors of desertification in Hunshandake Sandy Land from 1985 to 2017 
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年 国家 实施 的 农村 土地 承包 责任 制 和 1999 年 开始 
实施 的 退耕 还 林 还 草 、 围 封禁 牧 等 政策 相对 应 。 浑 
善 达 克 沙 地 人 口 呈 现 缓慢 平稳 增长 的 趋势 ,国内 生 
产 总 值 三 次 产业 变化 趋势 较 明 显 ,从 2005 年 开始 
出 现 明 显 增长 趋势 。 

气候 因素 对 沙漠 化 的 影响 是 一 个 漫长 的 过 程 ， 
因此 先 对 年 平均 气温 、 年 降水 量 、 年 平均 风速 和 年 
最 大 风速 等 自然 因子 进行 了 1~10 a 的 滑动 平均 值 
计算 ,并 与 DDI 之 间 进 行 了 关联 分 析 , 最 后 得 到 年 平 
均 气 温 、 年 降水 量 年 平均 风速 2a 累加 效应 和 年 最 
大 风速 8a 累加 效应 与 沙漠 化 程度 关联 性 较 高 。 
3.2.2 近 30a 浑 善 达 克 沙 地 沙漠 化 驱动 因子 分 析 
根据 11 个 驱动 因子 ,年 降水 量 (X1) 、 年 平均 气温 
(X2) 年 平均 风速 (Z3)、 年 最 大 风速 ( 双 )、 年 末 总 人 
O (X5) 耕地 面积 (X6) 载 冀 量 (X7)、 国 内 生产 总 值 
(X8)、 第 一 产业 (X9) .第 二 产业 (X10) 和 第 三 产业 
(X11) ,针对 浑 善 达 克 沙 地 近 30 a 的 因子 数据 进行 
主 成 分 分 析 ,选取 特征 根 大 于 1, 累计 贡献 率 大 于 
85.0% 的 主 成 分 因子 进行 主 成 分 载荷 分 析 ,11 个 影 
响 因 子 中 ,前 3 个 主 成 分 累计 贡献 率 达 到 89.9% ,各 
主 成 分 的 贡献 率 分 别 为 56.6% 17.8% .15.4%( 表 3)。 

在 第 一 主 成 分 中 ,变量 的 贡献 率 为 56.6% ,其 中 
年 最 大 风速 (X4) \ 年 末 总 人 口 (X5)、 国 内 生产 总 值 
(X8) .第 一 产业 (X9) .第 二 产业 (X10) 和 第 三 产业 
CI) 的 载荷 较 高 。 显 然 , 第 一 主 成 分 是 以 人 为 因 
素 为 主 , 其 中 年 末 总 人 口 .国内 生产 总 值 和 三 次 产 
业 等 因子 对 沙漠 化 的 发 展 过 程 起 到 正 向 作用 。 第 
二 主 成 分 中 ,年 降水 量 (X1) ,年 平均 风速 (X3) RA 
量 (X7) 等 3 个 因子 贡献 较 大 ,其 中 年 降水 量 的 贡献 
最 大 ,表明 该 主 成 分 以 自然 因子 为 主 ;第 三 主 成 分 
中 ,年 平均 气温 (X2) 年 平均 风速 (X3) 等 自然 因子 
和 耕地 面积 (X6) 载 畜 量 (X77) 等 人 为 因子 贡献 较 
大 , 且 自 然 因子 和 人 为 因子 对 该 主 成 分 的 贡献 率 约 
相当 。 综 合 来 看 ,在 1985 一 2017 年 间 浑 善 达 克 沙 地 
沙漠 化 过 程 受 人 为 因素 影响 较 大 (第 一 主 成 分 贡献 
率 最 高 ) ,自然 因素 也 起 到 一 定 辅助 作用 。 


4 讨论 


浑 善 达 克 沙 地 作为 中 国 北方 荒漠 化 最 严重 的 
地 区 之 一 ,已 被 列 入 全 国 重点 生态 环境 建设 地 区 。 
文章 利用 1985 一 2017 年 浑 善 达 克 沙 地 5 期 遥感 影 


表 3 1985 一 2017 年 浑 善 达 克 沙 地 沙漠 化 驱动 
因子 主 成 分 分 析 ( 载 荷 和 矩阵 ) 
Tab.3 Principal component analysis of driving factors 
(load matric) of desertification in Hunshandake Sandy 
Land from 1985 to 2017 


变量 PCI PC2 PC3 
年 降水 量 (X1) -0.125 0.851 0.079 
年 平均 气温 (X2) 0.547 -0.286 0.612 
年 平均 风速 (X3) -0.542 0.553 -0.458 
年 最 大 风速 (X4) -0.957 0.134 -0.018 
年 末 总 人 口 (X5) 0.968 -0.013 0.160 
耕地 面积 (X6) 0.319 0.293 0.803 
载 畜 量 (X7 ) -0.476 0.694 0.498 
国内 生产 总 值 (X8) 0.930 0.248 -0.232 
第 一 产业 (X9) 0.960 0.217 -0.125 
第 二 产业 (X10) 0.898 0.290 -0.261 
第 三 产业 (X11) 0.917 0.263 -0.211 
特征 根 6.226 1.959 1.691 
贡献 率 /% 56.600 17.800 15.400 
累计 贡献 率 /% 56.600 74.400 89.800 


注 :PC1 .PC2 .PC3 分 别 为 第 一 、 第 二 ,第 三 主 成 分 。 


像 数 据 ,通过 eCognition 9.0 和 ArcGIS 10.3 等 专业 软 
件 ,进行 沙漠 化 指数 的 提取 ,并 利用 主 成 分 分 析 法 
研究 沙漠 化 变化 影响 因子 ,明确 了 近 30 a BRIA 
沙 地 沙漠 化 时 空 演变 过 程 及 其 驱动 机 制 。 近 30 a 浑 
善 达 克 沙 地 沙漠 化 土地 面积 以 2002 年 为 转折 点 呈 
现 先 增长 后 下 降 的 趋势 ,其 中 重度 和 极度 沙漠 化 土 
地 面积 明显 减少 ,总 体 上 由 恶化 变 好 转 。 银 山 * 通 
过 MODIS 数 据 进 行 了 1989 一 2007 年 浑 善 达 克 沙 地 
荒漠 化 动态 研究 ,发现 以 2000 年 为 转折 点 出 现 先 恶 
化 后 改善 的 趋势 ,本 文 结论 进一步 验证 了 该 结 

而 浑 善 达 克 沙 地 沙漠 化 变化 主要 受 人 为 因素 影响 ， 
自然 因素 也 起 到 辅助 的 作用 。 随 着 人 口 数 量 的 不 
斯 增长 , 农 牧 业 活动 变 得 频繁 ,1980 年 以 来 ,中 国 实 
施 土 地 联 产 承包 责任 制 , 各 地 出 现 大 范围 开 若 、 尾 
种 现象 ,一 方面 提高 了 农业 生产 产量 ,提高 了 国内 
生产 总 值 ,但 男 一 方面 盲目 的 开垦 和 不 合理 的 土地 
利用 ,使 得 生态 脆弱 的 农 牧 交错 区 域 生态 环境 不 断 
恶化 。1988 年 开始 浑 善 达 克 沙 地 耕地 面积 急剧 增 
加 , 比 1987 年 约 增加 了 479 km ,粮食 种 植 面积 扩 
大 ,草场 被 占 为 耕地 的 同时 ,草场 载 畜 量 也 明显 增 
加 ,同期 约 增加 19.32x10' 只 ,因此 受 人 口 .国内 生产 
总 值 和 产业 发 展 等 人 为 因素 的 影响 ,使 得 其 沙漠 化 
面积 不 断 增加 , 沙漠 化 程度 不 断 加 重 。 到 了 2000 
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年 , 围 封禁 牧 和 退耕 还 林 还 草 等 政策 的 实施 以 及 
“三 北 防护 林 ” 生 态 工程 的 初步 成 效 , 浑 善 达 克 沙 地 
载 畜 量 和 耕地 面积 急剧 减少 ,土地 承载 力 下 降 ,使 
得 土地 沙漠 化 明显 减缓 ,成 为 沙漠 化 土地 由 恶化 变 
好 的 转折 点 。 显 然 , 人 为 因素 在 土地 沙漠 化 逆转 过 
程 中 起 着 决定 性 的 作用 。 王 涛 等 中 研究 科尔沁 地 
区 沙漠 化 过 程 驱 动因 素 中 ,也 得 出 自 2000 年 以 后 国 
家 实施 的 生态 环境 保护 政策 ,加 强 对 区 域 生态 环境 
的 保护 和 重建 力度 ,使 得 林地 和 草地 面积 增加 ,土地 
沙漠 化 面积 明显 减少 ,加 之 载 畜 量 的 降低 ,减轻 了 对 
草场 的 破坏 ,沙漠 化 呈现 逆转 趋势 。 岳 喜 元 等 所 进 
行 科 尔 沁 沙 地 沙漠 化 风险 评价 中 也 得 出 相同 的 结 
论 。 自 然 因素 中 ,1980 年 以 来 浑 善 达 克 沙 地 年 平均 
气温 变化 不 是 很 大 ,范围 在 1.5~4.9 °C ,平均 值 
3.5% ;年 降水 量 波动 变化 较 大 ,范围 约 在 200~400 
mm 之 间 ,平均 值 约 290 mm,2000 年 后 呈现 波动 增 
长 趋势 ;年 平均 风速 和 最 大 风速 有 所 降低 , 且 年 降 
水 量 相对 于 年 平均 气温 对 土地 沙漠 化 的 影响 较 大 ， 
年 降水 量 的 增加 和 风速 的 降低 有 助 于 植被 的 生长 
和 土地 沙漠 化 的 稳定 ,在 一 定 程 度 上 对 土地 沙漠 化 
起 到 缓解 作用 。 这 与 李 春 兰 等 ”进行 的 气候 波动 
下 浑 善 达 克 沙 地 荒漠 化 动态 变化 分 析 结 果 是 一 臻 
Wo Ton ES FE ESE 2000 一 2014 年 间 浑 善 达 
克 沙 地 植被 覆盖 度 变化 ,得 出 NDVI 年 际 变化 与 降 
水 量 的 关系 较 密 切 ( 尼 达 到 0.75) ,这 从 另 一 方面 说 
明了 降水 量 对 土地 沙漠 化 的 缓解 作用 。 

在 浑 善 达 克 沙 地 沙漠 化 驱动 因素 中 ,年 末 总 人 
口 、 国 内 生产 总 值 三 次 产业 、 载 畜 量 .耕地 面积 入 
是 沙漠 化 的 主要 影响 因素 ,而 年 降水 量 的 增加 和 年 
平均 风速 的 降低 在 一 定 程 度 上 对 土地 沙漠 化 起 到 
缓解 作用 ,其 他 因子 对 沙漠 化 土地 影响 较 小 。 而 在 
浑 善 达 克 沙 地 沙漠 化 的 治理 过 程 中 ,控制 牲畜 头 
数 ,限制 耕地 面积 的 继续 增长 ,实施 退耕 还 林 还 草 
和 禁 牧 、 休 牧 等 政策 ,并 通过 现代 化 技术 提高 农 牧 
业 产 量 , 保 证 经 济 增长 的 同时 ,合理 开展 治 沙 防 沙 
工程 ,进一步 改善 浑 善 达 克 沙 地 沙漠 化 程度 ,提高 
区 域 生 态 环境 质量 。 


5 结论 


=A 


本 文 以 浑 善 达 克 沙 地 为 研究 区 ,利用 5 期 影像 
数据 ,进行 时 空 特征 分 析 的 基础 上 进一步 完成 沙漠 


化 驱动 机 制 研究 。 结 果 如 下 : 

(1) 时 间 上 , 近 30 a, 浑 善 达 克 沙 地 沙漠 化 土地 
以 2002 年 为 转折 点 总 体 上 呈现 先 增 加 后 减少 的 趋 
势 ,其 中 重度 .极度 沙漠 化 面积 先 增加 后 减少 ,分 别 
由 占 沙 漠 化 面积 的 29.7% 、16.0% 减 少 到 22.6% 、 
6.3% , 轻 度 和 中 度 沙漠 化 土地 面积 都 有 所 增加 ,分 
别 由 20.9% .20.6% 增 加 到 32.1% 39.0% ,在 整个 沙 
漠 化 土地 中 所 占 比 例 较 大 。 

(2) 空间 上 , 浑 善 达 克 沙 地 沙漠 化 土地 呈现 由 
西向 东 逐 渐 减 轻 的 趋势 ,其 中 重度 沙漠 化 土地 主要 
分 布 于 西部 ,极度 沙漠 化 土地 从 西部 以 南 逐 渐 往 西 
北方 向 移动 ,中 度 沙漠 化 土地 主要 分 布 于 中 部 地 
区 , 轻 度 沙 漠 化 土地 主要 分 布 于 东部 地 区 , 非 沙 漠 
化 土地 分 布 于 东部 和 东北 地 区 , 占 总 土地 面积 比例 
由 8.1% 增 加 到 13.6%。 

(3) 浑 善 达 克 沙 地 沙漠 化 的 发 展 受 人 为 因素 影 
响 较 大 ,自然 因素 对 其 也 有 一 定 辅助 的 作用 。 而 人 
为 因素 中 年 来 总 人 口 .国内 生产 总 值 三 次 产业 对 
主 成 分 的 贡献 率 均 高 于 0.9 ,而 载 畜 量 和 耕地 面积 贡 
献 率 也 分 别 达 到 0.694 .0.803 ,对 沙漠 化 土地 的 影响 
较 显 著 ; 自然 因素 中 年 降水 量 的 增加 和 风速 的 降低 
对 土地 沙漠 化 也 起 到 一 定 程度 的 缓解 作用 ,年 降水 
量 和 年 最 大 风速 对 主 成 分 的 贡献 率 分 别 达到 0.851、 
0.957。 而 1980 年 实施 的 土地 联 产 承 包 责 任 制 和 
2000 年 以 后 国家 实施 的 一 系列 生态 环境 保护 政策 
对 浑 善 达 克 沙 地 沙漠 化 土地 由 增长 到 缓解 变化 起 
到 关键 的 作用 。 
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Abstract: The Hunshandake Sandy Land, as one of the most serious desertification areas in northern China, has 
been listed as a national key ecological environment construction area. A series of policies to control sand dust 
sources and restore the ecological environment in the Hunshandake Sandy Land has been implemented, such as 
returning farmland to forest and grassland and forbidding and restricting grazing, so as to curb the development 
momentum of desertification and improve the regional ecological environment. However, the ecological 
environment of the Hunshandake Sandy Land is fragile and, coupled with significant climate change, 
unreasonable land use patterns, excessive reclamation, grazing, and other phenomena, has not completely 
improved, resulting in the ecological environment of the area still facing problems. Therefore, it is of great 
significance to study the process of desertification and its influencing factors. The desertification index was 
calculated and the spatial-temporal evolution characteristics of desertification were studied in five remote sensing 
imagery periods of the Hunshandake Sandy Land from 1985 to 2017 using eCognition 9.0 and an object-oriented 
computer automatic classification method; further, the driving mechanism of desertification was studied using 
principal component analysis. In this paper, the research period of the Hunshandake Sandy Land desertification 
has been extended, the data have been updated (desertification data for 2017 were added), and the driving 
mechanism research has been synchronized with the desertification cycle. The results show that the 
desertification land area of the Hunshandake Sandy Land increased and then decreased over the last 30 years, 
with a turning point in 2002. The area of mild desertification first decreased and then increased (from 42.1% to 
20.9% and then to 32.1%), mainly in the east. Moderate desertification tended to increase gradually (from 20.6% 
to 39.0% ) and is distributed in the central region. However, severe desertification decreased from 30.3% to 
22.6% . Extreme desertification first increased and then decreased, increasing from 7.0% to 16.0% and then 
decreasing to 6.3%, and is distributed in the west and southwest and has gradually moved to the northwest. 
Overall, the spatial distribution of the desertification degree decreases from west to east. The development 
process of desertification in the Hunshandake Sandy Land has been greatly influenced by human factors, while 
natural factors play an auxiliary role. Of the human factors, the contribution rate of total population, GDP, and 
three industries to the principal component is higher than 0.9, while the contribution rates of the livestock 
carrying capacity and cultivated land area have reached 0.694 and 0.803, respectively, having a significant impact 
on desertification. Increasing precipitation and decreasing wind speed can alleviate desertification to some extent. 
The contribution rates of precipitation and the annual maximum wind speed to the principal component are 0.851 
and 0.957, respectively. The land contract responsibility system implemented in 1980 and a series of ecological 
and environmental protection policies implemented by the state after 2000 play a key role in the Hunshandake 
Sandy Land switching from growth to mitigation. Therefore, in the process of desertification control, we need to 
implement policies of returning farmland to forest and grass and grazing prohibition and restriction to decrease 
the desertification of the Hunshandake Sandy Land and to promote the healthy development of the regional 
ecological environment. 
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